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\ Abstract 

> ■ 

\f} ■ The Gutzwiller trace formula links the eigenvalues of the Schrodinger operator H 

as Planck's constant goes to zero (the semiclassical regime) with the closed orbits of 
the corresponding classical mechanical system. Gutzwiller gave a heuristic proof of this 
■ trace formula, using the Feynman integral representation for the propagator of H. Later, 

using the theory of Fourier integral operators, mathematicians gave rigorous proofs of 
the formula in various settings. Here we show how the use of coherent states allows us 

^ 1 to give a simple and direct proof. 

-S : 

1 Introduction 

0: 

Our goal in this paper is to give a simple proof of the "semiclassical Gutzwiller trace 
formula" . The pioneering works in quantum physics of Gutzwiller [|l7]] (1971) and Balian- 
Bloch || (1972-74) showed that the trace of a quantum observable A, localized in 
a spectral neighborhood of size 0{Ti) of an energy E for the quantum Hamiltonian H, 
can be expressed in terms of averages of the classical observable A associated with A 
over invariant sets for the flow of the classical Hamiltonian H associated with H. This 
is related to the spectral asymptotics for H in the semi-classical limit, and it can be 
understood as a "correspondence principle" between classical and quantum mechanics 
as Planck's constant h goes to zero. 

Between 1973 and 1975 several authors gave rigorous derivations of related results, 
generalizing the classical Poisson summation formula from d 2 /d9 2 on the circle to elliptic 
operators on compact manifolds: Colin de Verdiere ||, Chazarain Q, Duistermaat- 
Guillemin [14]. The first article is based on a parametrix construction for the associated 



heat equation, while the second two replace this with a parametrix, constructed as 
a Fourier integral operator, for the associated wave equation. More recently, papers 
by Guillemin-Uribe (1989), Paul-Uribe (1991, 1995), Meinrenken (1992) and Dozias 
(1994) have developed the necessary tools from microlocal analysis in a nonhomogeneous 
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(semiclassical) setting to deal with Schrodinger-type Hamiltonians. Extensions and 



simplifications of these methods have been given by Petkov-Popov [31|, and Charbonnel- 
Popov ||. 

The coherent states approach presented here seems particularly suitable when one 
wishes to compare the phase space quantum picture with the phase space classical flow. 
Furthermore, it avoids problems with caustics, and the Maslov indices appear naturally 
In short, it implies the Gutzwiller trace formula in a very simple and transparent way, 
without any use of the global theory of Fourier integral operators. In their place we use 
the coherent states approximation (gaussian beams) and the stationary phase theorem 

The use of gaussian wave packets is such a useful idea that one can trace it back 
to the very beginning of quantum mechanics, for instance, Schrodinger (1926). 
However, the realization that these approximations are universally applicable, and that 
they are valid for arbitrarily long times, has developed gradually. In the mathematical 
literature these approximations have never become textbook material, and this has lead 
to their repeated rediscovery with a variety of different names, e.g. coherent states 
and gaussian beams. The first place that we have found where they are used in some 
generality is Babich M (1968) (see also ||). Since then they have appeared, often as 
independent discoveries, in the work of Arnaud [||] (1973), Keller |^4| (1974), Heller [2C] 



(1975, 1987), Ralston @,[f4| (1976,1982), Hagedorn g§] (1980-85), and Littlejohn || 
(1986) - and probably many more that we have not found. Their use in trace formulas 
was proposed by Wilkinson |36| (1987). The propagation formulas of []l8| were extended 
in Combescure- Robert with a detailed estimate on the error both in time and in 
Planck's constant. The early application of these methods in [0] was for the construction 



of quasi- modes, and this has been pursued further in |33j| and Paul-Uribe |g7y. There 



have also been recent applications to the pointwise behaviour of semiclassical measures 



Acknowledgement: The authors thank J. Sjostrand for helpful discussions of this 
topic, and J. Ramanathan for valuable comments on the preliminary version of this 
paper. 

2 The Semiclassical Gutzwiller Trace Formula 

We consider a quantum system in L 2 (R n ) with Hamiltonian 

H = -h 2 A + V(x), (1) 

where A is the Laplacian in L 2 (R n ) and V{x) a real, C°°(]R n ) potential. The corre- 
sponding Hamiltonian for the classical motion is 

H(q,p)=p 2 + V(q), 
and for a given energy E (e H) we denote by Eg the "energy shell" 

^:={(W)£R 2 " : H(q,p)=E}. (2) 
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More generally we shall consider Hamiltonians H obained by the ft-Weyl quantization 
of the classical Hamiltonian H, so that H = Op£(H), where 

Ofl(H)rl,(x) = (2ich)^J H ^e^^dydi (3) 

The Hamiltonian H is assumed to be a smooth, real valued function of z = (x, £) G R 2n , 
and to satisfy the following global estimates 

• (H.O) there exist non- negative constants C, m, C 7 such that 

\d]H(z)\ < C 7 < H(z) >, Vz G R 2n , V 7 G N 2n (4) 
< fl-(^) > < C < H(z') > ■ < z - z' > m , Vz, z' G R 2n (5) 

where we have used the notation < u >= (1 + In) 2 ) 1 / 2 for u G M m . 

Remark 2.1 i) H(q,p) = p 2 + V(q) satisfies (H.O), ifV(q) is bounded below by some 
a > and satisfies the property (H.O) in the variable q. 

ii) The technical condition (H.O) implies in particular that H is essentially self-adjoint 
on L 2 (R n ) for h small enough and that x(H) is a h-pseudodifferential operator if x G 
C °°(R) (see ^) . 



Let us denote by <f>± the classical flow induced by Hamilton's equations with Hamil- 
tonian H , and by S(q,p;t) the classical action along the trajectory starting at (q,p) at 
time t = 0, and evolving during time t: 

S(q,p;t) = [ (p s ■ q s - H(q,p))ds (6) 
Jo 

where (qt,Pt) = 4>t(Q,p), and dot denotes the derivative with respect to time. We shall 
also use the notation: at = (j>t{cx) where a = (q,p) G R 2n , is a phase space point. 

An important role in what follows is played by the "linearized flow" around the 
classical trajectory, which is defined as follows. Let 



H"(a t ) (>2JI 



da 2 



(7) 



a=at 



be the Hessian of H at point at = <j>t(&) of the classical trajectory. Let J be the 
symplectic matrix 

H-i o) 

where and / are respectively the null and identity n x n matrices. Let F(t) be the 
In x In real symplectic matrix solution of the linear differential equation 

f F(t) = J H"{a t ) F(t) 

I \ W 
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F(t) depends on a = (q,p), the initial point for the classical trajectory, a t . 

Let 7 be a closed orbit on J2e with period T 7 , and let us denote simply by F y the 
matrix = F(T 7 ). F 7 is usually called the "monodromy matrix" of the closed orbit 
7. Of course, F 7 does depend on a, but its eigenvalues do not, since the monodromy 
matrix with a different initial point on 7 is conjugate to F 1 . F 7 has 1 as eigenvalue of 
algebraic multiplicity at least equal to 2. In all that follows, we shall use the following 
definition 

Definition 2.2 We say that 7 is a nondegenerate orbit if the eigenvalue 1 of F 1 has 
algebraic multiplicity 2. 

Let a denote the usual symplectic form on R 2n 

a(a, a') = p • q - p • q a = (q,p); a' = (q',p') (10) 

(• is usual scalar product in R"). We denote by {a^o^} the eigenspace of F 7 belonging 
to the eigenvalue 1, and by V its orthogonal complement in the sense of the symplectic 
form a 

V = {aeR 2n : a(a,a ± ) = a(a,a[) = o} . (11) 

Then, the restriction P 7 of F 1 to V is called the (linearized) "Poincare map" for 7. 

In some cases the Hamiltonian flow will contain manifolds of periodic orbits with the 
same energy. When this happens, the periodic orbits will necessarily be degenerate, but 
the techniques we use here can still apply. The precise hypothesis for this ( "Hypothesis 
C") will be given in Section 4. Following Duistermaat and Guillemin we call this a "clean 
intersection hypothesis" , but it is more explicit than other versions of this assumption. 
Since the statement of the trace formula is simpler and more informative when one does 
assume that the periodic orbits are nondegenerate, we will give that formula here. 

We shall now assume the following. Let (Te)t be the set of all periodic orbits on J2e 
with periods T 7 , < |T 7 | < T (including repetitions of primitive orbits and assigning 
negative periods to primitive orbits traced in the opposite sense). Then we require: 

• (H.l) There exists SE > such that H- l ([E -6E,E + 5E]) is a compact set of 
M 2 " and E is a noncritical value of H (i.e. H(z) = E => VH{z) / 0). 

• (H.2) For any T > 0, (r^T is a discrete set, with periods — T < T 7l < ••• < 

T 1N < T. 

• (H.3) All 7 in (Te)t are nondegenerate, i.e. 1 is not an eigenvalue for the corre- 
sponding "Poincare map", P 7 . 

We can now state the Gutzwiller trace formula. Let A = Op^(A) be a quantum 
observable, such that A satisfies the following 

• (H.4) there exists 5 G R, C 7 > (7 G N 2n ), such that 

\d]A(z)\ < C 7 < H(z) > S W G R 2 ", 

• (H.5) g a C°° function whose Fourier transform g is of compact support with 
Supp g C [-T,T] 
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and let x be a smooth function with a compact support contained in ]E — 5E, E + 
SE[, equal to 1 in a neighborhood of E. Then the following "regularized density 
of states" Pa(E) is well defined 



PA{E) = Tr (x(H)A X (H)g j j (12) 

Note that (H.l) implies that the spectrum of H is purely discrete in a neighborhood of 
E so that pa(E) is well defined. Then we have the following, 

Theorem 2.3 : Assume (H.0)-(H.3) are satified for H, (H.4) for A and (H.5) for g. 
Then the following asymptotic expansion holds true, modulo 0(h°°), 

PA (E) = (7r)- n / 2 g(0)^ {n ' 1] I A(a)da E (a) + £ c k (g)h k (13) 

J ' Ee k>~n+2 

7 e(r E ) T [ l aet ^ ^Jl Jo j>i 

where A{a) is the classical Weyl symbol of A, 

T* is the primitive period of '7, 

cr-y is the Maslov index of 7 ( 07 G 2 ), 

Sry = §^ pdq is the classical action along ~f, 

c k(d) are distributions in g with support in {0}, 

cfj(g) are distributions in g with support {T 7 } and da E is the Liouville measure on J2e : 
do~E= |^^| {dTiE is the Euclidean measure on Tie) 

Remark 2.4 We can include more general Hamiltonians depending explicitly in h, 

K 

H = WH® such that ffW satisfies (H.0) and for j > I, 
i=i 

\SPH^{z)\ < C 7d < HW(z) > (14) 

It is useful for applications to consider Hamiltonians like + fiH^ where 
may be, for example, a spin term. In that case the formula ( |l~3| ) is true with differ- 
ent coefficients. In particular the first term in the contribution of T 7 is multiplied by 

exp(-i/ T * HW(a s )ds). 



Remark 2.5 For Schrodinger operators we only need smoothness of the potential V . In 
this case the trace formula is still valid without any assumptions at infinity for V 
when we restrict ourselves to a compact energy surface, assuming E < limin^i^^ V(x). 
Using exponential decrease of the eigenf unctions Jill/ we can prove that, modulo an error 
term of order 0(h +0 °), the potential V can be replaced by a potential V satisfying the 
assumptions of the Remark (2.1). 
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3 Preparations for the Proof 

We shall make use of "coherent states" which can be defined as follows. Let 



M%) = W n/4 exp J-HL-J , (15) 

be the ground state of the n-dimensional harmonic oscillator, and for a = (q,p) 6 H 2n , 

T(a) = exp i^~(p ■ x - q ■ KD X )\ (16) 

is the Weyl-Heisenberg operator of translation by a in phase space where D x = -M^ We 
also denote by 

<p a = T(a)ipo (17) 

the usual coherent states centered at the point a. Then it is known that any operator 
B with a symbol decreasing sufficiently rapidly is in trace class (see fl5|]), and its trace 
can be computed by 

Tr£ = (2nhy n J < tp a , Bip a > da. (18) 
The regularized density of states pa(E) can now be rewritten as 

p A (E) = {2n)- n - 1 n- n J g{t) e iEt/h < ip a ,A x U(t) ip a > dtda (19) 

where U (t) is the quantum unitary group : 

U(t) = e - itS / h (20) 

and A x = x(H)A X (H). 

Our strategy for computing the behavior of pa(E) as % goes to zero is first to compute 
the bracket 

m(a, t) =< A x (p a , U(t)(p a >, (21) 
where we drop the subscript \ m A x for simplicity. It is useful to rewrite ( |i~5|) as 
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= AftV'o, (22) 
where is the following scaling operator: 

(AfcV) (x) = h~ n ^ ^ (xfi- 1 / 2 ) and fo(x) = vr^exp (-\x\ 2 /2) . (23) 

First of all we shall use the following lemma, giving the action of an ?i-pseudodifferential 
operator on a Gaussian. 

Lemma 3.1 Assume that A satisfies (H.0). Then we have 

Aip a =5>^" T)Q + 0(h°°) (24) 

7 ^' 



2n 2n 

in L 2 (R n ), where 7 e N 2r \ | 7 | = 7! = n 7j ! and 

1 1 

= T(a)A h O P y(z^ . (25) 

where Opf^z' 7 ) is the 1-Weyl quantization of the monomial : 
(x,^ = ^",7 = (V,7")e]N 2 «. 

This lemma is easily proved using a scaling argument and Taylor expansion for the 
symbol A around the point a. Thus m(t, a) is a linear combination of terms like 

m y (a,t) =< * 7>a , U(t)ip a > . (26) 

Now we compute U(t)<p a , using the semiclassical propagation of coherent states result as 



it was formulated in Combescure-Robert [ 10 1 . We recall that F(t) is a time dependent 



symplectic matrix (Jacobi matrix) defined by the linear equation (||). Met-F denotes the 
metaplectic representation of the linearized flow F (see for example Folland flj|), and 
the ^-dependent metaplectic representation is defined by 



Meta(F) = A" 1 Met(F)A n (27) 



We will also use the notation 



S(a, t) = r Ps • q s ds - tH(a) - Pt ' q \ P ' q (28) 
JO 2 

From Theorem (3.5) of |l(J (and its proof) we have the following propagation estimates 
in the L 2 -norm: for every N G IN and every T > there exists Cn,t such that 

\\U(t)<p a -exp T(a t )Meta( J P(t))A a Pjv(^,-Da ; ,t^)^|| < C^r^ (29) 



where P/v(i, ^) is the (h, t)-dependent differential operator defined by 

P N (x,D x ,t,h) = I+ J2 K k/2 - j Pkj(x,D,t) 

with I N = {(k,j) £NxN, 1 < j < 2N - 1, fe > 3j, 1 < fe - 2j < 2N} (30) 

where the differential operators Pkj(x, D x ,t) are products of j Weyl quantization of 
homogeous polynomials of degree k s with J2i< s <j = k (see [10] Theorem (3.5) and 
its proof). So that we get 

Pkj( x i D x,t)il)o = Qkj(x)i/> (x) (31) 

where Qkj(x) is a polynomial (with coefficients depending on (a,t)) of degree having 
the same parity as k. This is clear from the following facts: homogeneous polynomials 
have a definite parity, and Weyl quantization behaves well with respect to symmetries: 
Op w {A) commutes to the parity operator T,f(x) = f(—x) if and only if A is an even 
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symbol and anticommutes with S if and only if A is an odd symbol) and ifto(x) is an 
even function. So we get 

m{a, t) = 2^ cfej, 7 ft 2 J exp I — - — 

(j,k)eI N ;\j\<2N V 

T(a)A h Q~fc,T(a t )A n Q k Mtf,(F(t))fo) + O(^) (32) 



where Qk,j respectively Q 7 are polynomials in the x variable with the same parity 
as k respectively |7|. This remark will be useful in proving that we have only entire 
powers in fi in fll3|), even though half integer powers appear naturally in the asymptotic 
propagation of coherent states. By an easy computation we have 

T (a)A h Q^ , T (a t )A h Q k>j Met(F(t))^ 
exp ( - a t)) (Ti (^—7^) Q^o, Q fc ,iMet(F(i))Vo ) (33) 



where T\{-) is the Weyl translation operator with fi = 1. 
We set 

m kdl7 (a,t) = l^i(^^Q^,Q k jMet(F{t))^ (34) 

m (a,t) = ^7i(^^)^,Met(F(t)M^ (35) 

We compute mo (a, ) first. We shall use the fact that the metaplectic group transforms 
Gaussian wave packets to Gaussian wave packets in a very explicit way. If we denote 
by A, B, C, D the four n x n matrices of the block form of F(t), 

F{t) = (i *) (36) 



C D 

it is clear, since F is symplectic, that U = A + iB is invertible. So we can define 

M = VU- 1 , where V = (C + iD). (37) 

We have Ch.4) 

m (a,t) = (det C/)~ 1/2 ^ n/2 ■ 
ia n GXP + iI ) x • x ) - ~ 9 2 qt ) ■ (p-Pt + i(q - qt))j dx (38) 



Remark 3.2 In ( z (0)c the following meaning: if t \—> z(t) is a continuous 

1/2 

mapping from IR into C \ {0} such that z(0) > then (z(t)) c denotes the square root 
defined by continuity in t starting from y/z(0) > 0. Thus factor (det U) c 1 ^ 2 in (38) 
records the winding of det U(t) at t varies. This takes the place of the "Maslov line 
bundle" in this construction. 



<s 



If we make the change of variables x i-> (y — qt)j\ffi in (32) and hence in (38), then 
the formula for the regularized density of states in (19) takes the form 

p A (E)=[dt[ da [ a(t,a,y,h)e^ E< - y ' a ^dy. (39) 
Jr JR? n JR n 

The phase function &e is given by 

®E(t,y,a) = 

S(a, t) + q-p + {y - q t ) -Pt + \{y ~ qt) ■ M(t)(y - q t ) + % -\y - q\ 2 - y ■ p + Et, (40) 

where • denotes the usual bilinear product in C n , and a = (q,p), at = <pt(&) as before. 
Our plan is to prove Theorem 2.3 by expanding (39) by the method of stationary phase. 
The necessary stationary phase lemma for complex phase functions can easily be derived 
from Theorem 7.7.5 in [22, Vol. 1]. There is also an extended discussion of complex 
phase functions depending on parameters in [22] leading to Theorem 7.7.12, but the 
form of the stationary manifold here permits us to use the following 

Theorem 3.3 (stationary phase expansion) Let O C M d be an open set, and let 
a,f€ C°°(0) with 9/ > in O and supp aCO. We define 

M = {x€O,Xf(x) = 0,f'(x) = 0}, 

and assume that M is a smooth, compact and connected submanifold of IR d of dimension 
k such that for all x G M the Hessian, f"(x), of f is nondegenerate on the normal space 
N x to M at x. 

Under the conditions above, the integral J(to) = f R d e w ^ x ^a{x)dx has the following 
asymptotic expansion as uj — > +oo, modulo 0(uj~°°), 



d-k 

2n\ — 



CjU) 



-j 



j>0 



The coefficient cq is given by 



j 

JM 



det 



f"(m)\N m \]-^ 2 

a(m)dVM{m), 



where dVM(m) is the canonical Euclidean volume in M, mo £ M is arbitrary, and 

— 1/2 

[detP]* denotes the product of the reciprocals of square roots of the eigenvalues of P 

f" (rri)\N 

chosen with positive real parts. Note that, since Of/ > 0, the eigenvalues of — ^-p— - lie 
in the closed right half plane. 

Sketch of proof : Using a partition of unity, we can assume that O is small enough 
that we have normal, geodesic coordinates in a neighborhood of M. So we have a 
diffeomorphism 
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where U is an open neighborhood of (0, 0) in M k x M d h , such that 

X(x',x") eM^x" = 
and if m = x( x ', 0) G M we have 

X f (x>,0)(R k x )=T m M 

x'(x',0)(Rp k ) = N m M, normal space at m G M). 
So the change of variables x = x(x',x") gives the integral 

J( w ) = ^e i ^^( a; '' a; "))a(x / ,x")|detx / (x / ,x")|dx / dx // . (43) 

The phase 

/(*',*") := fix&y)) 

clearly satisfies 

{£„(x', x") = 0, 3/V, x") = 0} ^ x" = 0. (44) 

Hence, we can apply the stationary phase Theorem 7.7.5 of [22], (Vol. 1), in the variable 
x", to the integral (43), where x' is a parameter (the assumptions of [22] are satisfied, 
uniformly for x' close to 0). We remark that all the coefficients Cj of the expansion can 
be computed using the above local coordinates and Theorem 7.7.5. 

4 The stationary Phase Computation 

In this section we compute the stationary phase expansion of (39) with phase $e given 
by (40). Note that a(t,a,y,H) is actually, according to (32), a polynomial in h 1 ^ 2 and 
hT 1 ! 2 . Hence the stationary phase theorem (with fi independent symbol a) applies to 
each coefficient of this polynomial. 

The first order derivatives of <&E(i, y, a) (up to 0((y — q) 2 , (a — at) 2 ) terms) are given 

by 

dt $e = E - H(a) +(y- q t ) • Pt ~ Qt • M(y - q t ) 
d y <&e =Pt~P + i{y - q) + M(y - q t ) 

d q $ e = i(q - q t ) - l A{p-p t ) + ( *C - l AM - il){y - q t ) 
d p <S> E = q-q t + ( t D- f BM - I)(y - q t ). 

Furthermore, since F is symplectic, one has 

21m $e = \v- q\ 2 + \(A + %B)~\y - q t )\ 2 . 
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This implies that is critical on the set 



C E = {(y,a,t) £ IR n y x m 2 a n x M t : y = q t , a t = a, H{a) = E}. 
Thus each component .M 7 of Ce has the form 

M n = {(?/, a, t) = (q, a, T(a)) : a = (p,q) G 7, a T(a) = a, #(a) = e} . (45) 

We will assume that each 7 is a smooth compact manifold. One sees immediately that 
the manifolds 7 are unions of peridic classical trajectories of energy E. We will also 
assume a "clean intersection" hypothesis which we will state shortly. Thus we have 
assumed that 

C E = {0}xEU{M 71 ,...,M 1N }, (46) 

where each has the form (45) with 7^ in the fixed point set of the mapping a 1— ► ax k ■ 
The first thing to check, in order to apply the stationary phase theorem is that the 
support of q in (39) can be taken as compact, up to an error 0(h°°). We do this in 
the following way: let us recall some properties of ft-pseudodifferential calculus proved 



in [21, [L2|. The function m{z) =< H(z) > is a weight function. In |13] it is proved 



that x(H) = H x where H x G S(m k ), for every k. More precisely, we have in the H 
asymptotic sense in S(m~ k ), 

j>0 



and support [H X J\ is in a fixed compact set for every j (see (H.5) and [21| for the 
computations of H x ,j)- us recau that the symbol space S(m) is equipped with the 
family of semi-norms 

sup m (z)\— — u(z)\ 



Now we can prove the following lemma 

Lemma 4.1 There is a compact set K in IE? n such that for 

m(a, t) =< A x ip a , U (t)(p a > 

we have 

f \m(a,t)\da = 0(h + ™), 

JIR 2n /K 

uniformly in every bounded interval in t. 

Proof: Let x € C^(]E - 5E,E + 5E[) such that XX = X- Using (H.4) and the 
composition rule for ?i-pseudodifferential operators we can see that A X (H) is bounded 
on L 2 (IR n ). So there exists a C > such that 

\m(a,t)\<C\\ x(H)tp a f • 

But we can write 
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II x(H)ip a || 2 =< x(H) 2 cp a ,ip a > . 

Let us introduce the Wigner function, w a , for <p a (i.e. the Weyl symbol of the orthogonal 
projection on (p a ). We have 



< x(Hy<p a ,<p a >= (irh) n J H x 2{z)w a {z)dz 

where 

_ z-a| 2 

w a (z) = (ttH) n e & . 
Using remainder estimates from pl[ we have, for every N large enough, 



H X 2= H x 2 >h h j + h N+1 R N (h) 

0<j<N 

where the following estimate in Hilbert-Schmidt norm holds 

sup || Rn(K) \\hs< +00. 
o<h<i 

Now there is an R > such that for every j, we have Supp[i/ X 2 j] C {z, \z\ < R}. So 
the proof of the lemma follows from 

II Rn{K) \\hs= (2vr^)- n j || R N (h)ip a || 2 da, 
and from the elementary estimate, which holds for some C, c > 0, 

\a-r\ 2 _c 

e s dzda < Ce s . □ 

|«|<-R,|a|>-R+l 

The next step is the computation of the Hessian of $b on a ■M.^k- After an easy 
but tedious computation, with the variables ordered as (t,y,p,q), the Hessian <fr" E is the 
following (1 + 3n) x (1 + 3n) matrix 



cb" 



( H p ■ (JET, + MH p ) -H q - H P M -H p (D - MB) -H P (C - MA) \ 
—H q — MHp M + il D-MB-I C-MA-il 

-(*£>- t BM)H p t D- l BM-I l BMB - l DB l BM A — f BC 
\ -(*C- t AM)H p t C- l AM-iI l AMB — f CB l AMA- l CA + iI ) 

(47) 

where H p (resp. H q ) denotes the vector d p H\ a=at (resp. d q H\ a=at ), A, B, C, D, are the 
n x n matrices given by (36), 1 A the transpose of A, and M is defined by (37). (Recall 
/ is the identity matrix). 

We are going to perform elementary row and column operations on (47) to compute 
the nullspace of <& E , and the determinant of & E restricted to the normal space to the 
critical manifold. To begin with we have H\ = t R^' E RQ where 
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and Hi is given by 



( I \ 
H p I 
0/0 

V o oo/; 



Hi = 



( H p • (-H g + iH p ) 
-H q + iH p 



-Hg + iH p 



-H n 



V 



—Hp 
-iH„ 



M + il D- MB - I 

l D- l BM-l f BMB — l DB 
t C- l AM-iI t AMB — f CB 



-iHp \ 



C - MA - il 
l BMA — f BC 
f AMA — l CA + iI ) 



Multiplying H\ on the right by 



( 1 \ 

/ B A 

0/0 

\ I J 



changes it to 



/ H p (-H q + iHp) 



Ho 



V 



'Hg + iHp 

—Hp 
-iH n 



-H q + iHp 

M + il 
t D- f BM — I 
*C- f AM-iI 



-Hp + (-H q + iH p )B {-H q + iHp) A - iH p \ 
D-I + iB C + i(A-I) 

-B I- A 

-iB -i(A-I) 



The key simplification comes from (37) which gives M = (C + iD)(A + iB) 1 , and 
hence, since F is symplectic 

t D- l BM =[ t D(A + iB) - ^{C + iD)](A + iB)- 1 = (A + iB)- 1 

t C- l AM =[ t C(A + iB)- l A{C + iD)](A + iB)- 1 = -i(A + iB)-\ 

Thus, subtracting the appropriate multiples of the third row in Hi from the other rows 
we get 



/ 





-H q 

-Hr, 



\ 



—Hp 
D-I 
-B 




-H q \ 



C 
I- A 




(-Hg + iHp^A + iB)- 1 
(C + iD + U)(A + iB)- 1 
(/ - A - iB){A + iB)' 1 
^(A + iB)- 1 

Finally using the fourth row to remove the three upper entries in the second column, 
multiplying the third row by —1, interchanging the second and fourth rows, and the 
third and fourth columns, we arrive at the simple form 



H A 






Hp 
-H n 





^{A + iB)- 1 





— Hg —Hp \ 





A- I 

C 




B 

D-I J 



(48) 
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and H4 = i?i$gi?2 where R\ and R2 can be computed by repeating the elementary row 
and column operations that we have performed on the identity matrix, and in particular 
det Ri = 1 and det R 2 = (-l) n . 

In order to apply the stationary phase theorem the null space of <&" E must be the 
tangent space to the critical set Ce- However, one can read off the null space of H4 from 
(48) 



Null H 4 = i^Null 



(t,0,v,w) : (F-I)\ 



+ r 



H p 
-H n 



and H q ■ v + Hp ■ w = 



(49) 



This leads us to impose the following "clean flow condition" 



Hypothesis C: Assume that T>e '■= {(01, t) G x IR /<l>t{ot) = a} is a submanifold of 
IR l+2n . Then we say the T>e satisfies the clean flow condition, if for any (a,t) G T>e, 
the tangent space to T>e is given by 



T a , t V E = I (v, w, r) G IR 1+2n :(F-I) Q 



+ r 




= and H q ■ v + H p ■ w = o| . 

(50) 



Since Ce = {(y, a, t) : (a, t) G T>e and y = q}, the tangent space T y ,o,,t Ce equals 

H n 



{(t,v,w,v) :(F-I)\ 



~H n 



= and H q ■ v + Hp ■ w = 0}, 



and, assuming Hypothesis C, this does equal the null space of & E , since 

/ 

Av + Bw + tH p 



Ri = 



I * \ 





v 

V w J 



V 



w 

V 





( 


T 


\ 






V 








W 




) 


\ 


V 


I 



for (r, v,w) as in (49). Therefore, if P denotes the orthogonal projection on the null 
space of 3>£, then det($ E + P) will be the determinant of the Hessian of the phase 
restricted to the normal space, and setting 

P = R1PR2 (51) 

we have det(i^4 + P) = — (— l) n de,t(Q" E + P). Hence the computations of our paper 
provide a proof for the existence of a Gutzwiller trace formula under Hypothesis C. 
However, as stated earlier, we will only carry out the computations for the case that 7 
consists of a single trajectory here. In this case Hypothesis C reduces to the assumption 
(H.3) of isolated nondegenerate periodic orbits, and we may complete the computation 
in the following way. 

To compute det (H4 + P) we will use a special basis B. We denote by E\ the (alge- 
braic) eigenspace of F belonging to the eigenvalue A. Then under assumption (H.3) 
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dim © E x = 2n - 2 
dim E- 



l = ^ 

and a(E\,Ei) = for A / 1 where a is the symplectic form, as in (10). Let (21,22) be 
a basis for £1 with 

z 1 = {2H 2 p + H 2 q )- 1 / 2 {H p -H q ), 
and (F — I)z 2 = &z\. Let mi, • • • m2 n -2 be a (real) basis for the span of © E\, and let 
eo, • • • e„ be the Euclidean basis for R n+1 . Then we take B to be the basis 

{(e , 0) • • • (e n , 0)} U {(0, mi) ■ ■ ■ , (0, m 2ri _ 2 )} U {(0, (0, z 2 )}- 

Since the vector -P(^) spans the range of P and #4^) = 0, we can use column opera- 
tions to remove the contribution of P from all columns of the matrix H^ + P with respect 
to B, except the one corresponding to Z\. Then we can use column operations to remove 
all entries in the z\- and z 2 -columns corresponding to the basis vectors (ei, 0) • • • (e n , 0), 
and (0, mi) • • • (0, m 2 „_ 2 ). Note that this does not change the entries in the first row of 
the matrix, since a(zi,mj) = 0, j = 1, . . . 2n — 2. After these simplifications which do 
not change the determinant, the matrix of H4 + P with respect to B becomes: 



( 








b \ 





-2i(A + iB)' 1 














Pj-I 





V a 








n J 



The vector a is just ((2H 2 + H 2 ) 1 / 2 , 0, • • • 0) and 

b= (x,-(2H 2 + H 2 ) 1 / 2 a(z 1 ,z 2 
Therefore the determinant of H4 + P equals 



det 



A + iB 



-1 



det(P 7 -7) detO, 



where 



(52) 



(53) 




x 
(2^ + F 2 )V2 x 

c 



-(2H 2 + H 2 )a( Zl ,z 2 ) 




(54) 



Here x is used for entries that do not enter the calculation, and c is the component of 
P(0, z\) along the basis vector z 2 . 
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To compute c and finish the computation of the determinant, we first compute 
P(0, z\). Writing z\ = (v,w), we have 



1 ^ 




( 


X 


\ 





1 




X 




V 


~ 2 




—w + Bv — Aw 




V w J 




K 


3v — Cw + Dv 


I 



(55) 



We let P\Z\ denote the last 2n components of P(0, z{). Since l FJz\ = Jz±, the normal- 
ization in the definition of z\ gives, cr(zi, P\z{) = 1. Therefore, if Pi^i = CZ2 + d^i we 
clearly have c = o"(zi, z-i)~ x . Thus (54) yields 



det^ = -(2H 2 + H 2 ) 
and, combining this with (53) and (56), we have 



(56) 



U 



det^u (A , 7) = (-ir- i Hrdet V2 

Using (42) and (57), we conclude 



\(0,H p ,-H q ,H p )\ 2 -det(P 7 -I). (57) 



4 = 



iS-y/H 



T* 



-l^KO^^^det^-J) 



det (|) 



-1/2 



det- 



CA-i/2 



A(a fl )dF(s). 



Using |(0, p s , g s ,p s )| _1 dU(s) = ds we get the result for dg in (13). Since det("P 7 — I) = 
(— l) a | det("P 7 — where a 1 is the number of real eigenvalues of V~, which are greater 
than 1, we see that 



fdet^-J =±i n - l+a '\det(T> y -I)\- l/2 . (58) 

* 

Note that the role of the Maslov index in (13) is to determine the sign in (58) and <r 7 
in (13) is either n — 1 + a' or n + 1 + a'. 

The other coefficients, dj are spectral invariants which have been studied by Guillemin 
and Zelditch. In principle we can compute them using this explicit approach. This com- 
pletes the proof of Theorem 2.3. 



(-l) 1 ~"det(P 7 -I) 
~ det (2) 
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